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Abstract 

We develop a closed-form three-parameter model for unsaturated hydraulic conductivity associated with 
a three-parameter lognormal model of moisture retention, which is based on lognormal grainsize distribu- 
tion. The derivation of the model is made possible by a slight modification to the theory of Mualem. We 
extend the three-parameter lognormal distribution to a four-parameter model that also truncates the pore 
size distribution at a minimum pore radius. We then develop the corresponding four-parameter model for 
moisture retention and the associated closed-form expression for unsaturated hydraulic conductivity. The 
four-parameter model is fitted to experimental data, similar to the models of Kosugi and van Genuchten. 
The proposed four-parameter model retains the physical basis of Kosugi's model, while improving fit to 
observed data especially when simultaneously fitting pressure-saturation and pressure-conductivity data. 
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1. Introduction 



Kosugil (|1994I ) assumed pore size is a lognormal random variable and derived a three-parameter model 



for moisture retention, the three parameters being the mean and variance of the pore-size distribution and 
the maximum pore radius. In the limiting case where the maximum pore radius a pproaches infinity the 



three-parameter model simplifies to a two-parameter model for which 



form expression for unsaturated hydraulic conductivity using the theory of 



Kosugil (11996 ) devel o ped the closed 



Mualeml (jl976| ). kosugil (|1996T ) 



did not develop a three-parameter closed-form equation for hydraulic conductiv ity, but reverte d to the two- 



parameter form, owing to diffic ulty in analytically integrating the expression of 



the work of 



Kosugi (1994 



Mualeml (|1976l ). We extend 



1996) and develop closed-form expressions for unsaturated hydraulic conductivity 



associated with the three-parameter lognormal moisture retention model. The derivation of the closed-form 
equation for un saturated hydraulic conductivity is made possible by a slight modificatio n to the theor y of 
Mualeml (|1976l ). Further, we modify the pore-size probability density function (PDF) of iKosugil (| 1994T ) by 
incorporating a nonzero minimum pore radius, as suggested by iBrutsaerti ([1966). This modification results 
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in a four-parameter moisture retention model and a corresponding four-paramet er closed-form e quation for 
unsaturated hydraulic conductivity, again obtained using the modified theory of iMualeml (|1976l ). The four 
paramet ers in th e prop osed model are all based on physical properties of the porous medium (similar to the 



model of 



Kosugil (|1996l) ). not fitting parameters without physical significance. 



2. Theory 



The lognormal distr ibuti on is commonly use d to statistically characterize pore size in granular porous 



media. 



BrutsaertJ (|1966l ) a ndlKosue 



(1994 



19961 ) considered lower- and upper-tail truncated lognormal PDFs 



for pore-size distributions. iBrutsaerti (|1966f ) considered the log-transformed random pore radius i?— ro, where 



ro is the radius at which the effective moisture content vanishes (associated with residual saturation). 
2.1. The three-parameter lognormal model 



The classical (non-truncated) lognormal distribution for pore size (jBrutsaertl . 119661) . measured here by the 



random pore radius R £ [0, oo]. For a real p orous medium R £ [0, r max ], where r max is some finite maximum 



pore radius. To correct for the finite interval. 



Kosugil (|1994T) used the random variable R e = (1/R — 1 /r max )~ 



to rescale the classical PDF. The PDF of R is related to the PDF of R e according to 



MO 



Sr. 



(1/r - l/r max ) 1 



(1) 



(1 - r/r max j 

According to Young-Laplace theory, the capillary pressure head h and pore radius r are related according 
to h = k/t, where k = 2*y cos a/ (pg), 7 is interface surface tension, a is the interface contact angle, p is fluid 
density, and g is gravitational acceleration. For water in a glass tube, k « 0.149cm 2 . It then follows that 
the PDF of the random capillary pressure head H is 



/ile (h/n-l/r max y 
[h/K - l/r max ) 

Further, if R e is lognormally distributed, the PDF for the capillary pressure head H can be written as 



Mh) 



1 



■ exp 



log (h - h c ) - n v 



V2< 



(Jz 



(2) 



(3) 



<2no- z (h-h e ) 

for all h > h c , where h c — n/r max is the bubbling pressure he ad, = log(At ) — \iz is the mean of log(iJ), o 2 z 
is the variance of fiz is the mean of Z, and Z — log (R e ). iKosugil (|1994f ) used the dimensionless random 
variable R' e = -R e max to obta in the PDF in Q with the parameters p,z and o~ 2 z scaled appropriately. 
As shown by 



Kosugil (|l994l ). the three-parameter moisture retention curve that follows from (|3|) is 

io g (h-M-M h > hci 

I h < h r . 



'(h) 



lerfc ( - 



(4) 



where 9*(h) = (9(h) — 9 r )/(9 s — 6 r ) is moisture capacity, 9(h) is volumetric moisture content, 9 r is residual 



42 moistu re content, 9 S is saturated moisture content, and erfc is the complementary error function. 



Kosugi 



(1994) did not dev elop a corresponding closed- form equation for unsaturated hydraulic conductivity. 



Mualeml (|1976l) developed a functional relation between unsaturated hydraulic conductivity K(9*) = 



K s K r and capillary pressure head h 



K, 



da; 
~hjx) 



o h(x) 



(5) 



where K r and K s are relative and saturated hydraulic conductivity and x is an integration variable. Equation 
([5]) can be rewritten in terms of the PDF of capillary pressure head as 

-i 2 



K r 



1 



d.r 



f H (x) 



dx 



(6) 



Using the theory of I Mualeml (|l976[ ) , iKosugil (jl996l) develo ped the two-p arameter closed- form equation for 



unsaturated hydraulic condu ctivity by forcing r 



oo. 



52 the theory of 



Kosugil (|l996l ) made this simplification because 



Mualeml (|l976l ) as given in (JS]) is not readily amenable to integration when r max in the three- 



parameter lognormal distribution is finite. 

We obtain a clos ed-form express ion for unsaturated hydraulic conductivity for finite values of 
modifying ([5]) due to iMualeml (|l976l ) into 

2 



Kr 



d.r 



h(x) 



d.r 



h(x) 



by 



(7) 



based on the assumption fn(h)/h « fn(h) /(h— h c ). Using this approximation the closed-form unsaturated 
hydraulic conductivity expression for the three-parameter lognormal model is 

2 



K r (h) 



■{» 



erfc 



lQg(fc-fec)-Mt?+Cz 



1 



h > h c 
h < h r 



(8) 



In the limit as r n 



derived by 



oo, ([8]) reduces to the two-parameter unsaturated hydraulic conductivity expression 



Kosugil (|l996l ). Figure [T] shows that the truncated lognormal pore-size distribution ([S]) results in 



the conductivity curves shifted to the right (K r = 1 at h = h c ) compared to the curves corresponding to the 



two-parameter model of 



Kosugil (|l996l ). where K r = 1 is reached at h > 0. 



64 2.2. Four-parameter lognormal model 



To incorporate the lower-tail truncation of iBrutsaertl (|l966l ) into the distribution of IKosugil (|l994T I. we 
introduce the random variable R e defined by 

4=(V —) , (9) 



R-r 



3 



which can be shown to yield the following lognormal PDF for capillary pressure head, 



Mh) = 



1 



2iru<Jz 



exp 



log(u) - fi. n 
V2< 



(10) 



for all h £ [h c , ft m ax] where u — (1/h — l/h max ) — h c and h max = is the pressure head associated with 
the smallest undrainable pores in the medium corresponding to the physical boundary on pore radius. Figure 
[5] shows the three lognormal PDFs for capillary pressure head: the classical (non-truncated) lognormal dis- 
tribution, the upper-truncated lognormal dis tributio n and the doubly truncated lognormal distribution 



10]) . It can be seen in the figure the model of 



Kosugi (1994) departs from the classical lognormal distribution 



only at small head values (large pore radii) whereas the proposed four-parameter distribution departs from 
the classical distribution at both the lower and upper limbs of the function. 



The moisture retention curve is derived from (|10p in a manner similar to that done by iKosugil (|19961 ) for 
the four-parameter lognormal distribution, and is given by 



9*(h) 



icrfc 



lQg(M)-^,, 



h c < h < h n 
h < h c , 



(ii) 



Finally, it can be shown that the closed-form expression for unsa t urate d hydraulic conductivity using the 
doubly truncated PDF (fTD|) and the modified equation of . 



'{ierfc 



log(M)-^,,-Q- 



Mualeml (|l9761 ) © is 

2 



K r (h) 
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h ^ ^max- 



(12) 



In the limit as both r„ 



from 



Kosugi (1994, 



1996) 



oo and r — > 0, (ITT1) and (|T2^) reduce to corresponding two-parameter expressions 



3. Fitting four-parameter lognormal model to experimental data 

The four- parameter logno r mal m odel for moi s ture r etention (fTTj) was fitted to experimental data (same 



Kosugil (|199q) ) using non-linear least squares to estimate the pa- 



van Genuchtenl (|1980l ) and 

/iz and a\. We also estimated 9 r and 9 S from experimental data, rather than using 



87 data used by 

sa rameters r , r max , 

89 reported values. Using estimated parameters in the four-parameter model (|12|) . predictions of unsaturated 

90 hydraulic conductivity were compared to measured values. The results for Hygiene sandstone and silt loam 

91 G.E. 3 are shown in Figures [3] and 2] A summary of estimated parameters for these soils is given in TablcQ] 

92 Figure EJb) shows a solid curve computed using the approximate, but closed- form expression for K r ([7]) and 

93 a dashed curve numerically integrated from the traditional Mualem relationship ([S]), whi ch are nearly iden- 

94 tical for this set of parameters. The fits to experimental data are comparable to those of lKosugil (|199q) and 
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van Genuchtenl (|198(J ). Markov-chain Monte Carlo simulations with the four-parameter model and moisture 
retention curve data show that the logarithms of parameters r$ and r max are sometimes identifiable only in 
a threshold sense. No well-defined optimal value exists, the data imply the log minimum pore size is less 
than a threshold, or the log maximum is greater than a threshold. 

The four-parameter lognormal model was also fitted to moisture retention data for Beit Netofa clay, 
using these estimated parameter values to make predictions of conductivity. The four-pa rameter mode l 
yield ed predictions of uns aturated conductivity comparable to the two-parameter models of iKosugil |1996) 



and 



van Genuchtenl <| 19801 ) as shown by the dotted curves in Figure^a) and (b) (labeled 2- and 3-parameter). 



However, when hydraulic conductivity and moisture retention data were jointly used to estimate the param- 
eters ro, r maj; , [xz, c^, r , and 9 S , there was marked improvement in model fits to the data as shown by the 
solid curves in Figure EJa) and (b) (labeled 4-parameter). The parameter values used in the fitted model 
for moisture retention data shown in Figure |5ja) are the same values used for the fitted conductivity model 
in Figure [S^b). The s ame procedure did not yield similar results when applied to the the two-parameter 
model of lKosugil ((1996), as shown by the dashed curves in plots (a) and (b) of Figure which were obtained 



by jointly using moisture retention and conductivity data to estimating the parameters fiz and <7 2 Zl with 



— > oo. 



4. Discussion 

The proposed four-para meter lognormal mode l fits moisture ret ention data from three representative 



soils similar to the models of 



van Genuch tcn (1980) and 



conductivity are comparable to those of 



van Genuchten 



1980) and 



Kosugj (| 1996 ). Predictio ns of unsaturated hydraulic 



Kosugil (11996) for Hygiene sandstone and 



silt loam G.E. 3. Estimating model parameters for Beit Netofa clay using only moisture rete ntion data does 
not yi eld g ood prediction s of unsaturated hydraulic conductivity (similar to the models of 



van Genuchten 



(|l980l ) and IKosugil (|199q) ). It was essential to use both moisture retention and hydrauli c conductivity data 



us to estimate the param eters and improve the fit of the proposed model over the models of 



van Genuchten 



Kosugil (|l994j ) and 



(1980). This indicates moisture retention data alone may not be sufficient to estimate all 
four (or six, when including 9 r and 9 S ) model parameters. If conductivity measurements are available, they 
should be used with moisture retention data, to arrive at a more realistic closed-form model for unsaturated 
soil moisture retention and hydraulic conductivity. 

One would expect a model with a larger number of adjustable parameters to fit observed data better than 
similar models with less parameters, but an improved fit is often at the expense of the physical significance of 
the parameters. The proposed four-parameter model does no t make its gain s in model-data fit at the expense 



Kosugil <| 19941 ) used in deriving his model: they 



of parameter realism. The parameters follow the philosophy 
are related to pore-size di stributio n stati stics, rather than being fitting parameters. The development and use 
of the modified model of iMualeml (|1976f ) enables one to obtain approximate but closed-form expressions for 
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129 
130 
131 



unsaturated hydraulic conductivity for both the three- and four-param eter lognormal pore size distributions. 
The model derived here is a generalization of the the lognormal model of lKosugil ([19961 ) for non-zero minimum 
pore radius and non-infinite maximum pore radius. 
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Table 1: Estimated parameter values for proposed 4- parameter model 



Medium 


r (m) 




fmax(m) 




crz 


9 S 


r 


Hygiene Sandstone 


1.07 x 10 


-4 


2.52 x 10" a 


-6.300 


0.337 


0.250 


0.153 


Silt Loam G.E. 3 


1.48 x 10 


-4 


1.27 x 10~ 2 


-7.927 


1.118 


0.395 


0.192 


Beit Netofa Clay 


4.36 x 10 


-7 


1.32 x 10" 2 


-11.06 


2.333 


0.450 


0.100 
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classical 

proposed 

Kosugi (1994) 




h (cm) 

Figure 2: Comparison of lognormal capillary pressure head PDFs for the classical, three-parameter {Kosugi, 1994), and proposed 
four-parameter d 1 1 ) models. 




Hygiene sandstone 
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Figure 3: (a) Moisture content data and fitted moisture re tention curve, a nd (b) measured and predicted K r (using both 
analytical and numerical integration) for Hygiene Sandstone jMualeml <1976ft soil index 4130). 



8 




0.5 



0.4 



0.3 



0.2 



- 1 1 - 




Beit Netofa clay 



data o 

2- parameter 

3- parameter - 

4- parameter 



J L_ 



10° 10 1 10 2 10 3 10 4 10 5 

h (cm) 



10 u 



10" 



- 10 



10" 



-2 



10" 



10 



-5 



1-1 — I 1 — 



-©- 




(b) 



data o 

2- parameter — 

3- parameter - 

4- parameter 



10° 10 1 10 2 10 3 10 4 10 5 

h (cm) 



Figure 5: (a) M oisture content data and fitted moisture retention curve, and (b) measured and fitted K r for Beit Netofa Clay 
(Mualcm (1976)) soil index 1006). 
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